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§ ■ ABSTRACT 
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D ■ We present a quantitative estimate of the confusion of cluster radio halos 
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and galaxies in the measurement of the angular power spectrum of the thermal 

Sunyaev-Zel'dovich (SZ) effect. To achieve the goal, we use a purely analytic 

approach to both radio sources and dark matter of clusters by incorporating 

^ ■ empirical models and observational facts together with some theoretical consid- 

^ ■ erations. It is shown that the correction of cluster radio halos and galaxies to the 

CN ■ 

. measurement of the thermal SZ angular power spectrum is no more than 20% at 

/ > 2000 for observing frequencies z/ > 30 GHz. This eliminates the concern that 
. the SZ measurement may be seriously contaminated by the existence of cluster 

^ . radio sources. 

Q ■ Subject headings: cosmology: theory — diffuse radiation — radio continuum: 

general 

> 

^ ■ 1. Introduction 

As a result of gravitationally driven shocks and compression, majority of the baryons 
in the most massive, collapsed objects in the universe such as groups and clusters of galax- 
ies exist in the form of hot plasma with temperature of 10'' — 10^ K. These hot baryons 
manifest themselves by strong diffuse X-ray sources in the sky within the framework of 
bremsstrahlung emission, which are directly detectable with current X-ray instruments. On 
the other hand, the hot electrons in the virialized massive dark halos also scatter the passing 
cosmic microwave background (CMB) photons, giving rise to a subtle change in the CMB 
spectrum, known as the Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 1972). Unlike 
the X-ray emission of the hot gas in massive halos, which shows a strong dependence on the 
detailed structures and enrichness history of the hot gas as well as the distances of the host 
halos from us, the SZ effect is completely determined by the total thermal energy of the gas 
intrinsic to the systems. Therefore, it is well suited for studies of distant, massive clusters. 
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which may provide a stringent constraint on cosmological parameters especially for the mean 
matter density Qu and the normalization of matter power spectrum (Tg (e.g. Barbosa et al. 
1996; Molnar, Birkinshaw & Mushotzky 2002; Komatsu & Seljak 2002). A more practical 
and powerful utilization of the SZ effect nowadays is perhaps to statistically measure the 
angular power spectrum of the SZ signals at various scales. This allows us to gain the weak 
SZ signals at high statistical significance levels without resolving individual groups and clus- 
ters. Indeed, recent detection of the excess power relative to primordial CMB anisotropy 
at arcminute scales has been attributed to the statistical signals of the thermal SZ effect 
(Dawson et al. 2001; Bond et al. 2003). A number of ongoing and future experiments will 
soon be able to provide a more robust measurement of the SZ angular power spectrum at 
high multipoles or small angular scales. 

There are many other effects that should be corrected for before the SZ angular power 
spectrum is used for cosmological purposes. These include nongravitational effect on the 
distribution and evolution of the hot gas inside halos (da Silva et al. 2000, 2001; Seljak, 
BurwcII & Pen 2001; Xue & Wu 2001; Holder & Carlstrom 2001; Zhang, Pen & Wang, 2002; 
Komatsu & Seljak 2002; Zhang & Wu 2003; etc.), possible asymmetric shape of dark halos 
(Cooray 2000; Lee & Suto 2003; Piffaretti, Jetzer & Schindler 2003), contaminations of radio 
point sources (Cooray et al. 1998; Lin, Chiueh & Wu 2002; Holder 2002; Rubino-Martin & 
Sunyaev 2002), etc. Despite large uncertainties, many attempts have been made over the 
past few years towards the evaluation of these effects on the SZ cluster counts and angular 
power spectrum measurements. In this work, we will explore the contamination of radio 
sources in clusters in the measurement of the thermal SZ angular power spectrum, namely, 
cluster radio halos and galaxies. 

In theory of hierarchical clustering, massive dark halos like clusters of galaxies form by 
gravitational aggregation of individual low-mass objects. Sub-halo mergers will generate a 
considerably large number of energetic electrons that can be accelerated by cluster magnetic 
fields over cluster regions, giving rise to the radio halos due to synchrotron emission (Buote 
2001 and reference therein). Regardless of the debate on the acceleration mechanism for 
relativistic electrons in clusters, about 1/3 of nearby, rich clusters exhibit radio halos, and 
this fraction may increase with redshift if mergers were more frequent in the past (EnBlin & 
Rottgering 2002). Furthermore, there is a strong observational evidence that cluster radio 
halos show a morphology similar to the diffuse X-ray emission from the hot intracluster 
gas, and the corresponding radio power is correlated strongly with the X-ray luminosity 
and temperature of clusters (Liang et al. 2000; Govoni et al. 2001a,b). This enables one to 
estimate the distribution and evolution of cluster radio halos based on the well-determined X- 
ray properties of clusters. For example, one can easily derive the radio luminosity function of 
clusters at different wavebands and redshifts (EnBlin & Rottgering 2002). Together with the 
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radio surface brightness profiles of clusters motivated by empirical formula (e.g. (5 model) for 
the X-ray surface brightness distribution of clusters, one can estimate the radio background 
from the diffuse radio halos of clusters. 

Existence of diffuse radio halos constitutes a major source of contaminations in the 
measurement of the thermal SZ effect toward clusters. Since radio halos always contribute a 
positive signal superimposed on the CMB sky, the decrements of the CMB signals due to the 
thermal SZ effect of the hot intracluster gas at low frequencies would be compensated, while 
the increments of the SZ signals at high frequencies will be correspondingly enhanced. The 
crossover point occurs &i v = 217 GHz. Moreover, the amplitude of such a modification to 
the thermal SZ signals depends critically on observing frequencies. This arises because the 
radio flux of cluster halos as a result of synchrotron emission has a steep spectrum of roughly 
Sn oc Consequently, at very high frequencies the cluster radio halos may become hardly 
visible, and no modification to the SZ signals from the cluster radio halos is needed. However, 
with the decrease of observing frequencies the SZ decrements and the diffuse radio emission 
of cluster halos may appear to be comparable, and eventually the latter would become the 
dominant component in the microwave sky. This will lead to a nonnegligible correction to 
the SZ measurement (EnBlin 2002) especially in the central regions of clusters, because the 
radio emission is more centrally concentrated than the thermal SZ effect if the radio emission 
is assumed to follow the X-ray emission. Therefore, it deserves a close investigation of the 
infiuence of the cluster radio halos on the measurement of the SZ angular power spectrum. 

Another source of contaminations in the SZ measurements is radio point sources, which 
include radio galaxies inside clusters and tracing large-scale structures of the universe. In 
a similar way to cluster radio halos, cluster radio galaxies also compensate the decrements 
of the thermal SZ intensity in the Rayleigh- Jeans region. On the other hand, the Poisson 
distribution and clustering of radio galaxies on large scales would enhance the radio fiuctua- 
tions, resulting in a positive contribution to both primary CMB and SZ power spectra. The 
role of these field radio galaxies has been well addressed by several authors (e.g. Pranceschini 
et al. 1989; Toffolatti et al. 1998), and the resulting contamination in the measurement of 
CMB anisotropics is found to be generally small or can be somehow removed. Nevertheless, 
a significant fraction of the SZ signals could be diluted by the unsubtracted radio sources 
inside clusters (Cooray et al. 1998; Lin, Chiueh & Wu 2002). Using a toy model of ra- 
dio point source subtraction, in which a certain central region is removed from all clusters. 
Holder (2002) showed that the true SZ power at I > 1000 can be underestimated by as much 
as 30% at = 30 GHz. A more reahstic model for the radial distribution of radio galaxies 
in clusters is thus needed to refine this estimation. 

In this paper, using a purely analytical approach, we perform a quantitative study of the 
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contamination of cluster radio halos and galaxies in the measurement of the angular power 
spectrum of the thermal SZ effect. It is hoped that such a study may clarify the issues of 
to what extent one has to take the correction of cluster radio sources into account in the 
cosmological applications of the thermal SZ angular power spectrum. Throughout the paper 
we will work with a flat cosmological model of ilu = 0.35 and Q\ = 0.65. 



2. Thermal SZ angular power spectrum 

The change in CMB temperature Tcmb in direction due to the thermal SZ effect of 
the hot electrons in massive clusters is simply 

= 9{x)y{e); (1) 



y{0) = J rieCTT ^ 



dx; (2) 



9{x) = ^^-4, (3) 

where x — hpu/kTcuB, and Te are the number density and temperature of electrons in 
clusters, respectively, and the integral is performed along the line of sight x- The angu- 
lar power spectrum of the temperature fluctuation AT{0)/Tcmb can be separated into the 
Poisson term C[ and the clustering term Cf: 

Ce = C[ + Cf. (4) 

Since the clustering term C^" is several orders of magnitude smaller than the Poisson term 
Cf at small angular scales I > 1000 where the SZ power becomes dominant in the CMB sky, 
we will neglect the contribution of source clustering to the SZ angular power spectrum. In 
the flat sky approximation, the Poisson term can be calculated by (e.g. Cole & Kaiser 1988) 

in which (PV/dzdft denotes the comoving volume per unit redshift and per steradian, and 
ye{M,z) is the Fourier transform ol y{9): 

ye^2nj y{e)Jo{ee)ede, (6) 

where Jo is the Bessel function of zero order. We adopt the standard Press & Schechter 
(1974) formahsm for the mass function of clusters 

d'N ^ _ [2 p 5,{z) da{M) f S'^jz) \ 
dMdV MnMa^iM) dM ^"""^ \ 2a\M))' ^ ' 
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where p is the mean mass density of the universe, 5c{z) is the hnear over-density of spherical 
collapse that virialized at redshift and (t{M) is the present rms mass fluctuation within 
the top-hat filter M = ATipR^/Z: 

a\M) = — / k^P{k)W\kR)dk, (8) 

in which W{x) — 3(sinx — xcosx)/x^ is the Fourier representation of the window func- 
tion. We parameterize the power spectrum of fluctuation P{k) oc k""'T'^{k), and take the 
fit given by Bardeen et al. (1986) for the transfer function of adiabatic CDM model T{k) 
with the shape parameter F = 0.21. The primordial power spectrum is assumed to be the 
Harrison-Zerdovich case = 1. The amplitude in the power spectrum is fixed using the 
rms fluctuation on an 8 h^^ Mpc scale, cxg, which will be taken to be cxg = 0.9 in this study. 



3. Contribution of cluster radio halos 

3.1. Properties of cluster radio halos 

Over the past few years, diffuse radio halos have been detected in a few tens of nearby, 
rich clusters. They often extend to a distance of ~ 1 Mpc from cluster centers, and have 
regular shape, low surface brightness and steep radio spectrum S^, oc Moreover, ob- 

servations suggest that there is a tight correlation between the radio power Pi.4ghz and the 
X-ray luminosity {Lx) of host clusters (Liang et al. 2000; Govoni et al. 2001b), which can 
be well described by the following fitting formula (EnBlin & Rottgering 2002) 

P1.4GHZ = a.lO- Watt Hz-^ {jo^^J^^) ' 

in which — 2.78, h,, — 1.94, and h^Q is the present Hubble constant in units of 50 km s~^ 
Mpc~^. In order to set up a hnk between the X-ray luminosity and total mass of the cluster, 
we first employ the best-fit Lx-T relation of Xue & Wu (2000), based on the hitherto largest 
sample of nearby X-ray clusters, 

\ in-0.032 ' '^-^ 



It has been shown that such an Lx-T relation exhibits no significant evolution at least out 
to 2; ~ 0.8 (Rosati, Borgani & Norman 2002 and references therein). We then convert the 
X-ray temperature into the total cluster mass using virial theorem (Bryan & Norman 1998): 

^^(^^Y" (11) 
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in which 

kT* = 1.39 keV fT[h^A,(z)E\z)Y/\ (12) 

h is the Hubble constant in units of 100 km Mpc~^, /t is the normalization factor, which 
will be taken to be /t = 0.8 in this work, and /^d^) is the overdensity of the virialized dark 
halos with respect to the critical density of the universe at redshift z. Equivalently, one may 
also employ the empirical Lx-M relation (Reiprich & Bohringer 1999) incorporating with a 
simple model of cosmic evolution for dark halos. Consequently, for a given halo of mass M 
one is able to estimate its radio power at 1.4 GHz through 

/L7-*\ 5.41 / ^ N 3.61 

P1.4GH. = 3.17 X 10^° h,,' Watt Hz-^ j^^j [Tq^J " ^^^^ 

Radio observations indicate that cluster radio halos usually have a steep power index, which 
is also well motivated by synchrotron emission process. We adopt a spectral power index of 
a = 0.8 (e.g. Giovannini et al. 1993) to convert the radio power at 1.4 GHz to the radio 
emission at other wavelengths. 

Furthermore, several observations have also reveled a similarity of the radio and X-ray 
morphologies in clusters (Deiss et al. 1997; Liang et al. 2000; Govoni ct al. 2001a; Feretti 
et al. 2001). This enables us to model the radio surface brightness of clusters using the 
empirical fitting formula for X-ray surface brightness profile, i.e., the /3 model, 

/ 2 \ -3/3+0.5 

A/RH(r) = A/rho (^1 + j ' (14) 

in which the normalization factor can be determined by the total radio power Px at a given 
wavelength A. Meanwhile, the electron number density in clusters follows 

/ ^2 \ -3/3/2 

rieir) = rieo ( ^ + ) ■ (1^) 

Here we introduce the universal baryon fraction f], = f2b/^M to fix the central electron 
density rieo- fh = Afgas/M, where Mgas is the total gas mass within rvk- The gaseous halo 
(also the radio halo) is truncated at the virial radius Tyir defined by 

M ^ —rl,Ac{z)pdz), (16) 

where Pc{z) is the critical mass density of the universe at z. We adopt a constant value of 2/3 
for the /? parameter, and assume that the core radius Tc of gas distribution is proportional 
to the size of the host halo: Tc oc Tyir- We specify the proportionality constant such that 
Tc = 0.15 Mpc for a massive cluster of M = IO^^H-^Mq at z = (cf. Komatsu & 
Kitayama 1999). 
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Finally, radio halos are preferentially found in clusters that demonstrate violent mergers 
and substructures, which provide an accelerating mechanism for the relativistic particles re- 
sponsible for radio emission (e.g. Buote 2001). Yet, the fraction of clusters that possess giant 
radio halos is still uncertain. Following Enfilin & Rottgering (2002), we adopt a constant 
radio halo fraction of /rh = 1/3. Nevertheless, we will also demonstrate the result for a 
fraction of unity for comparison. Recall that cluster mergers may occur more frequently in 
the past than at the present. 



3.2. Modification to the SZ angular power spectrum 

In the Rayleigh- Jeans region, the presence of radio halos can dilute the thermal SZ 
signals towards clusters. The significance of the effect depends critically on working fre- 
quencies. Recall that radio emission of clusters has a steep power spectrum of a ~ 0.8. At 
very low frequencies, the decrements of the CMB temperature due to the SZ effect could 
be completely compensated (EnBlin 2002) by cluster radio halos. Here we demonstrate the 
radio halo contamination using an observing frequency of = 30 GHz. The change in the 
CMB temperature is related to the intensity fluctuations by 

AT _ (e- - If AI ^^^^ 



where Jq = 2(A;Tcmb)^/(^pc)^. A/ is composed of two components, the thermal SZ signals 
A/gz and the cluster radio surface brightness A/rh, 

A7 = A7sz + A7rh, (18) 
or equivalently, we can separated the temperature fluctuation into 

AT ATsz , ATrh 



rj. rj. +7f7^- (19) 

-L CMB -L CMB -L CMB 

If the fraction of cluster radio halos is /rh, then the observed SZ power spectrum at small 
angular scales is 

j "^'^ I "^^^-ImT [(l-/RH)|^(x)y,(M,^)|2 + /RH|yK^,^)r] , (20) 

where y'i{M, z) is the Fourier transform of (ATsz + ATrh)/Tcmb- 

In order to highlight the effect of cluster radio halos on the SZ angular power spectrum, 
we demonstrate in Figure 1 the relative correction to the thermal SZ angular power spectrum 
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jAQl/Q, along with the contributions of clusters in different redshift ranges. At 1000 < 
I < 8000, the radio halo contamination in the SZ angular power spectrum arises mainly from 
clusters at 0.1 < z < 0.5. High- redshift clusters come into effect only at very small angular 
scales i > 10^, while large angular scales are dominated by low-redshift clusters. In Figure 2 
wc show the contributions of cluster radio halos with different cluster masses to the thermal 
SZ angular power spectrum. At large angular scales below i ~ 500, the radio contamination 
is governed by nearby, massive clusters, which can be easily removed from the SZ map because 
of their rarity. At smaller angular scales beyond i — 1000, the contamination arises mainly 
from moderately rich clusters at higher redshifts z > 0.2. Contributions from poor clusters 
and groups of masses less than 2 x 10^^ Mq are neghgible because of the absence of radio 
halos in these systems. To provide a quantitative illustration of the modifications due to 
cluster radio halos, we list in Table 1 the relative corrections |AC/|/Q for a set of i evaluated 
by Equation (20), assuming the radio halo fractions of /rh = 1/3 and /rh = 1) respectively. 
We have also shown in Table 1 how our estimates are affected by the uncertainties associated 
with the measurements of radio power. X-ray luminosities and total masses of clusters. To 
do this, we simply vary the radio power P^, X-ray luminosity Lx and normalization factor fx 
in the virial theorem by 20%, 20% and 30%, respectively, which are roughly representative 
of the uncertainties in current determinations of these quantities. It turns out that the 
modifications to the thermal SZ angTilar power spectrum from cluster radio halos are only 
minor, and the maximum correction is no more than 13% for £ > 2000 even in the extreme 
case of where /rh = 1 and /t- = 1. Major uncertainty in our calculation arises from the 
mass estimate of clusters characterized by the normalization factor fx, which could differ by 
a factor of ~ 2 among different simulations and observations (e.g. Bryan & Norman 1998; 
Huterer & White 2002). 



Table 1: Relative corrections to the SZ angular power spectrum from cluster radio halos 



/rh 


e = 2000 


e = 


3000 


i = 4000 


e = 5000 


1/3 


1 1 c+0.09+0.11+3.10(y* 
-'-•■'-^-0.14-0.29-1.03 




1-0.11+2. 18(v 
-0.23-0.81 


n «Q+0.07+0.11+1.48ci/ 
"-'•"^-0.11-0.21-0.70 


n c 0+0.07+0.12+0.930/ 
^•"-•"-0. 10-0.20-0.63 


1 


q 77+0.26+0.34+8.640/ 
•J-' ' -0.43-0.88-3.09 


'^•^"^-0.36- 


^0. 32+5. 78 (v 
-0.71-2.45 


9 /if;+0.22+0.34+3.64ci/ 
^•^"-0.32-0.64-2.11 


9 9f.+0. 22+0. 36+2.010/ 
^•^■^-0. 30-0. 59-1. 89 



* Listed in order are uncertainties arising from 20%, 20% and 30% variations in radio power. 



X-ray luminosity and mass estimate of clusters, respectively. 
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4. Contribution of cluster radio galEtxies 

4.1. Distribution of cluster radio galaxies 

A critical issue in the estimate of confusion of cluster radio galaxies in the measurement 
of the SZ angular power spectrum is to properly handle the distribution of cluster radio 
galaxies. This includes the luminosity function and radial profile of radio galaxies in clusters. 
We adopt the univariate radio luminosity function given by Ledlow & Owen (1996), based 
on a statistical complete sample of 188 radio galaxies in Abell clusters within z = 0.09 and 
at = 1.4 GHz. The univariate radio luminosity function /i.4ghz can be well fitted by a 
broken power law 

1 f 1 C?^radio/c^logPl.4GHz , li d ^01^ 

log /l.4GHz = log „ ^ = a + 6 log Pi.4GHz, (21) 

iVgal(0.3rAbcu) 

where A''gai(0.3rAbcii) is the total number of galaxies brighter than R = —20.5 within 0.3 Abell 
radius of the cluster center, and (a, b) are proportionality constants with (—0.150, 1.766) and 
(-1.433,33.667) for Pi.4ghz < lO^^'^ WKz^^ and Pi.4ghz > 10^'^'*^ WHz-\ respectively In 
order to convert the above radio luminosity function at = 1.4 GHz into the ones (/j,) at 
other frequencies, we use the spectral index distribution fitted by Lin, Chiueh, & Wu (2002) 
for a sample of 64 radio sources in 56 clusters observed by Cooray et al. (1998), in which 
the mean spectral index is a = 0.71. This allows us to take the dispersion of spectral index 
into account. The luminosity- weighted cluster radio luminosity function takes the form of 

^^^(P f\-! C^i(^^) +7i(^^)logP.; logP < logPb 

log [KU) - I ^^^^^ ^ ^^^^^ ^^g^^. ^^^p ^ ^^^p^ , [^Z) 

where 

Ci(i/) = 2.68 - 0.956-°°^^'^ 

71 (i/) = 0.85e-°-°°2^'' + 1.49 x lO'^u 

C2{v) = 30.42 + 3.296-°°^^'^ 

72(jy) = _o.43e-° °°'2'^ - 6.14 x 10"^, (23) 
and the break power Pb has a weak dependence on frequency: logPb(z^) = 23.86 — 0.00396i/. 

We assume that the number density profile of galaxies in clusters follows the King model 

/ 2\-3/2 

ngai(r)=ngai(0)(^l + ^j . (24) 

The surface number density is thus Sgai(r) = J n^a\{r)dx- A good approximation of the core 
radius is Tc — O.lryir- In order to fix the central density ngai(O), we introduce the so-called halo 
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occupation distribution, Nga.i{M) — J ngai(r)47rr^dr. We adopt the best fit analytic formulae 
of Shcth & Diaferio (2001) for the spiral [iVgai,s(^^)] and elliptical [A^gai,E(-^)] galaxies based 
on the GIF simulations (Kauffmann et al. 1999): 

A^gai,s(M) = (M/MsT' + o.5e-^['°s(^/^°"''^®)l' 
A^gai,E(M) = (M/ME)"^e-('^^'^"^®/^)' 

7Vg,i(M) = iVg,i,s(M) + iVg,i,E(M), (25) 

where Mg = 7 x W^h-^M^, as = 0.9, Me = 3 x W^H-^Mq, and qe = 0.75. 

Another parameter we need to fix before we proceed to our numerical calculation is 
A''gai(0.3rAbcii)) the galaxies brighter than R = —20.5 within 0.3 Abcll radius rAbcU- Instead 
of evaluating Ai"gai(0.3rAbeii) directly, wc turn to the fraction of galaxies with absolute mag- 
nitude brighter than -Rumit = —20.5, /gai. This can be achieved by employing the luminosity 
functions of galaxies in clusters, $5 and $e. We adopt the best fit Schechter forms in the 
r* band for early and late type galaxies derived from the Sloan Digital Sky Survey (Goto et 
al. 2002). The best fit parameters in Goto et al. (2002) are based on galaxies within 0.75 
h^} Mpc, which is only slightly larger than O.SrAbeii- We convert the r* magnitude to the 
standard Johnson-Morgan-Cousins system using r* = i?-|-0.16(y — it!)-|-0.13 (Fukugita et al. 
1996) and then assign the colors to early and late type galaxies according to the tabulated 
values given by Fukugita, Shimasaku & Ichikawa (1995). Moreover, we choose a lower mag- 
nitude cutoff of i?cut = —18.5 to guarantee the validity of the above adopted halo occupation 
distribution equation (25) which accounts for galaxies brighter than V = —17.7 -|- blogh. 
Taking all these factors into account, we have 



•^'"2— '^^ ^ •^'•e, limit ^ ^ 



gal 



^s{r*)dr* + ^E{r*)dr* ' 



(26) 



The total radio power of all galaxies in a cluster can be obtained by extrapolating the 
radio luminosity function of Ledlow & Owen (1996) determined with 0.3rAbeii to the whole 
cluster region out to virial radius, which yields 

A^gal(M) 



tot,;/ 



_iVgal(0.3r Abcll 
= /galA^gal(M) J PJMogPu. 



Pu{dN,^iJd log P,)d log 



(27) 



The total surface brightness of all radio galaxies measured at frequency u in direction 6 
toward a cluster of mass M and redshift z is thus 



A/rg(0) 



[1 + z)Ptot,u 

^^Dl{z) , 
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Sgal(^)/gal/P./.^logP. 

47r(l + 2;)3 



(28) 



where Da and Dl are the angular diameter distance and luminosity distance to the cluster, 
respectively. This adds to the SZ signals such that the observed radio intensity is A7 = 
A7sz + A/rg- The corresponding angular power spectrum reads 



where y'f'^M.z) is the Fourier transform of (ATsz + ATrg)/Tcmb, and the temperature 
fluctuation ATrg is related to A/rg following equation (17). 



In a similar way to Figure 1, we plot in Figure 3 the relative corrections to the SZ 
angular power spectrum \ACi\/Ci because of the presence of radio galaxies in clusters, 
together with the contributions of cluster radio galaxies in different redshift ranges. At 
/ > 1000, the distant radio galaxies a.t z > 0.2, which could be regarded as the unresolved 
radio sources, dominate the correction term. We list in Table 2 the relative corrections for 
a set of small scales, together with the uncertainties corresponding to 20% variation of the 
radio luminosity function. It turns out that the confusion of unresolved radio galaxies in 
clusters in the measurement of the thermal SZ effect is less than 7.5% on all scales beyond 
/ > 2000. This is somewhat smaller than the estimate (30%) of Holder (2002), who evaluated 
the effect of radio point sources by removing a certain central region (~ 1') from all clusters. 
The latter could simultaneously lead to an underestimate of the thermal SZ effect in the 
centers of all clusters especially for high-redshift clusters. 



Table 2: Relative corrections to the SZ angular power spectrum from cluster radio galaxies 




(29) 



4.2. Modification to the SZ angular power spectrum 



redshift ^ = 2000 ^ = 3000 ^ = 4000 ^ = 5000 



0.2<^< 5 4.39t[j|^% 4.63l°|o% 4.69t°|°% 4.69t°:i% 



o<.<5 6.1711:?;% 5.801?:^^% 5.iil2:^i% 5.i6lS:«I% 



5. Conclusions 



Contamination of radio sources has been a major concern in the measurement of thermal 
SZ effect, which include radio halos and galaxies of clusters, and clustering of radio galaxies 
tracing large-scale structures of the universe. The latter comes into effect only on large 
angular scales and appears to be relatively small as compared with the thermal SZ effect 
(Toffolatti et al. 1998). We have thus concentrated in this paper on the role of cluster radio 
halos and galaxies. The effect of cluster radio halos on the measurement of the SZ angular 
power spectrum has previously remained unknown, while the confusion of the cluster radio 
galaxies has been recently shown to be very significant (Holder 2002). Using an analytic 
model for both radio sources and dark matter of clusters, we have performed a quantitative 
analysis of the extent to which cluster radio halos and galaxies may contaminate the mea- 
surement of thermal SZ angular power spectrum. We have found that while the cluster radio 
sources can indeed cause some confusion in the measurement of SZ angular power spectrum, 
the relative correction |AQ|/Cj is well within 20% at I > 2000 for an observing frequency of 
u — 30 GHz. The modification becomes even smaller for high-frequency observations. We 
thus conclude that there is probably no need to take the confusion of cluster radio sources 
into account in the thermal SZ measurement, if our observing frequencies are chosen to be 
z/ > 30 GHz, unless the thermal SZ power spectrum is used for the purpose of precision 
cosmology. 
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Fig. 1. — Relative corrections of cluster radio halos to the thermal SZ angular power spectrum 
evaluated at = 30 GHz and for /rh = 1/3 (sohd line). Contribution of clusters in different 
redshift ranges are also shown. 
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